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ABSTRACT 

We report on spectroscopy of the white-dwarf companion of the millisecond radio 
pulsar PSR J1012+5307. We find strong Balmer absorption lines, as would be expected 
for a cool DA white dwarf. The profiles are much narrower than usual, however, and 
lines are seen up to H12, indicating that the companion has a low gravity and hence a 
low mass. This is consistent with the expectation — based on evolutionary considerations 
and on the mass function — that it is a low-mass white dwarf with a helium core. By 
comparing the spectra to model atmospheres, we derive an effective temperature T c g = 
8550 ± 25 K and a surface gravity \ogg = 6.75 ± 0.07 (cgs units). Using the Hamada- 
Salpeter mass-radius relation for helium white dwarfs, with an approximate correction 
for finite-temperature effects, we infer a mass Mwd = 0.16 ± 0.02 M Q . This is the lowest 
mass among all spectroscopically identified white dwarfs. We determine radial velocities 
from our spectra, and find a radial- velocity amplitude of 280 ± 15kms _1 . With the 
pulsar's radial-velocity amplitude, the mass ratio Mpsr/Mwd — 13.3 ± 0.7. From all 
constraints, we find that with 95% confidence 1.5 < M PSR /M < 3.2. 

Subject headings: binaries: close — stars: neutron — pulsars: individual (PSR J1012+5307) 
— white dwarfs 
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1. Introduction 

Knowledge of the properties of the white-dwarf 
companions of radio pulsars can provide unique con- 
straints on the characteristics and evolution of these 
binaries, as well as on those of its constituents. Of 
particular interest are mass determinations for mil- 
lisecond pulsars. Evolutionary models indicate that 
millisecond pulsars have acc reted up to 0.7 M Q (e.g., 
Van den Heuvel & Bitzaraki |l995| ). Hence, if the neu- 
tron star started with the "canonical" 1.4 Mq (e.g., 



Timmes, Woosley, & Weaver 1996 ; for a rece nt cen - 
sus, Van Kerkwijk, Van Paradijs, & Zuiderwijk 1995), 
it would now be > 2M Q . If so, this would strongly 
favor a stiff equation of state (EOS) at supranuclear 
density (e.g., Cook, Shapiro, & Teukolsky 1994). For 
softer EOS, like the one recently proposed by Brown 



& Bethe (1994), such a massive neutron star would 
collapse into a black hole. So far, constraints on the 
masses could only be set from Shapiro delay of the 
pulsed signal. As this requires (rare) fortuitous align- 
ment, a reasonable constraint is available for one sys- 
tem only, PSR B1855+09, for which Kaspi, Taylor, & 



Ryba p94|) found M PSR = 1.50j£^ M Q . 



The temperature of the white dwarf — combined 
with the mass — can be used to determine the age 
from theoretical cooling tracks. The cooling age gives 
a limit to the time since mass transfer to the neu- 
tron star ceased. With mass estimates from the mass 
function and approximate cooling models, such lim- 
its have been used to show that magnetic fields in 
neutron s tars d o not decay on a ~ 10 7 yr time scale 
(Kulkarni 1986| ), as had been argued previously (for 
a revie w, P hinney & Kulkarni 1994 ). Recently Bell 
et al. ( 1995 ) used a lower limit to the cooling age of 
the white-dwarf companion of PSR J0034— 0534 to in- 
fer that the pulsar's initial spin period was <J 0.6 ms. 
If so, this would also constrain the structure of neu- 
tron stars, excluding the stiffcr EOS and thus m aking 
it hard to have masses ;> 2 M Q (Cook et al. 1994 ; 
Haensel, Salgado, & Bonazzola 1995| ). Without ac- 
curate white-dwarf masses and good cooling models, 
however, such arguments remain inconclusive. 

Masses and temperatures of DA white dwarfs can 
be determined accurately from their spectra (e.g., 
Bergeron, Saffer, & Liebert 1992| ), and therefore we 
have started a spectroscopic investigation of all bright 
companions (see Van Kerkwijk |1996| ; currently, six 
companions with V <J 23 are known). First results 



Here, we describe spectroscopy of the white-dwarf 
companion of the millisecond pulsar PSR J1012+5307, 
discovered recently by Nicastro et al. (1995). The sys- 
tem's orbital period is 14.5 h, and the mass function 
indicates a low companion mass: Mwd = 0.11 Mq 
for Mpsr = 1.4 M Q and inclination i = 90°. At the 
pulsar's position, Nicastro et al. found a faint star on 
the Palomar sky survey. Lorimer et al. (1995) con- 



firmed the positional coincidence, and found that the 
star had magnitudes and colors consistent with those 
of a low-mass white dwarf at the distance implied by 
the pulsar dispersion measure (~0.5kpc). 

2. Observations 

Spectra of the companion of PSR J1012+5307 
(V = 19.5 mag) were taken at the Keck 10 m tele- 
scope with the Low- Resolution Imaging Spectrometer 
(LRIS). On 4 June 1995, the 300 line mm" 1 grating 
was used for a classification spectrum, covering the 
wavelength range of 3430 to 8460 A at 2.5Apix~ 1 . 
With a slit width of 1", the resolution was - 12 A. 
The conditions were photometric, and for calibration 
an exposure was taken through a 4" slit. Directly fol- 
lowing the object, the spectrophotometric standard 
Feige 34 (sdO, V = 11.23, Massey et al. |l988|) was 



observed through both slits, at very similar airmass. 

On 22-25 November 1995 and 21 January 1996, the 
600 line mm -1 grating and 0.7" slit were used to in- 
crease the resolution to ~4 A. The wavelength range 
covered was 3500-6000 A. During these runs, the con- 
ditions were not photometric, but Feige 34 was again 
observed, to obtain relative spectrophotometry and 
to verify the stability of the velocity determinations. 

The reduction of all spectra was done using MI- 
DASQ and programs running in the MIDAS environ- 
ment. The frames were bias-corrected, flat-fielded 
and sky-subtracted using standard procedures, and 
the spectra extracted using optimal weighting. Fluxes 
were determined by first calibrating the wide-slit ex- 
posure relative to that of Feige 34 (correcting for small 
differences in airmass using the extinction curve of 
Beland, Boulade, & Davidge 1988), and then scaling 
the narrow-slit exposure such that it had the same 
continuum flux. 



5 The Munich Image Data Analysis System is developed and 
maintained by the European Southern Observatory. 



can be found in Van Kerkwijk & Kulkarni (1995) 
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Fig. 1. — The spectrum of the companion of PSR J1012+5307. Shown in the left-hand panel are the spectrum 
taken with the 300 line mm -1 grating (top curve) and the average of the spectra taken with the 600 line mm -1 
grating (bottom curve; the spectra were shifted to zero velocity before averaging). The two were normalized to the 
observed V-band flux (V — 19.58; Lorimer et al. 1995), and are shown offset by +15 and —15 /iJy, respectively. Also 
shown are the observed broad-band fluxes and the best-fit pure Hydrogen model spectrum. The latter was derived 
from a fit to the profiles of H/3 up to H8, and has T c s — 8550 ± 25 K and logg = 6.75 ± 0.07. In the right-hand 
panel, the observed line profiles, including those of Ha and H9, are shown with the modeled ones superposed. 



3. Effective temperature, surface gravity and 
mass 

The spectra of the companion of PSR J1012+5307 
are shown in Fig. [j]. Strong Balmer lines are seen, 
from Ha up to H12. In field white dwarfs, the lines 
are usually broader, and the higher members of the 
Balmer series are rarely seen, as the corresponding 
excitation levels are quenched by the high pressure. 
This indicates that the white dwarf has a low surface 
gravity and hence a low mass. 

The line profiles could be reproduced well with a 
white-dwarf model atmosphere (Fig. [j]), with an ef- 
fective temperature T c g — 8550 ± 25 K and a sur- 
face gravity logg = 6.75 ± 0.07 (cgs units). ( For de - 
tails about the procedure, see Bergeron et al. 1992 .) 
By way of verification, we also fitted the broad- 
band fluxes determined by Lorimer et al. ( 1995 ) to 
a model atmosphere, with the surface gravity fixed to 
logg = 6.75. We find T cS = 8730 ± 110 K, consistent 
with the spectroscopic determination. 

The mass of the white dwarf can be estimated from 
the surface gravity using a mass-radius relation appro- 
priate for a white dwarf with a helium core (given the 
low mass, the star will never have burnt helium). Un- 
fortunately, helium white dwarfs have not been stud- 



ied as well as the usual carbon-oxygen white dwarfs, 
and only the zero-temperature mass-radius relation 
of Hamada & Salpeter (1961) is available (Table [j]). 
At 8500 K, however, a white dwarf with a logg as low 
as that for our white dwarf, will not have reached 
its zero-temperature configuration. To estimate the 



relative "puffiness", we have used the Wood (1995) 
models for white dwarfs with pure carbon cores and 
thick hydrogen atmospheres (Table Q). Inter- and 
extra-polating in the corrected mass-radius relation, 
our final estimate for the mass of the white dwarf is 
log M WD /M = -0.78 ± 0.06 (where the quoted un- 
certainty combines in quadrature an uncertainty of 
0.03 due to the uncertainty in logg, and an estimated 
0.05 due to the uncertainty in the mass-radius rela- 
tion). Hence, M WD = 0.16 ± 0.02 M@. 

4. Radial Velocities and the Mass of the Neu- 
tron Star 



We determined radial velocities for both the com- 
panion of PSR J1012+5307 and Feige 34 by fitting the 
H/3, H7 and US lines simultaneously to Lorentzian 
profiles, with the widths and velocities of the three 
lines forced to be the same. By way of verification, 
we also made fits with the velocities left free. The 
results were consistent with each other for all spectra 
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Table 1 

Mass-radius relation including 
finite- temperature effects 



M 


Ro 


log go 


-Rs50o/-Ro 


log g 


(M Q ) 


(Rq) 


(cm s~ 2 ) 




(cms -2 ) 


0.154 


0.0218 


6.95 


1.351 


6.69 


0.213 


0.0195 


7.19 


1.236 


7.00 


0.305 


0.0171 


7.46 


1.111 


7.36 


0.399 


0.0152 


7.67 


1.071 


7.62 


0.499 


0.0137 


7.86 


1.048 


7.82 


0.609 


0.0122 


8.05 


1.032 


8.02 


0.734 


0.0108 


8.24 


1.022 


8.22 


0.885 


0.0092 


8.46 


1.010 


8.45 



Note. — The first two columns give the zero- 
temperature mass-radius relation of Hamada & 
Salpeter (1961) for a helium white dwarf. The fourth 
column lists the adopted ratio of the radius at T = 
8500 K and the zero-temperature radius, as deter- 
mined from the cooling curves of Wood (1995) for 
carbon white dwarfs with thick hydrogen envelopes. 



Table 2 

Radial Velocity Measurements 



PSR .71 01 2-1 


-5307 


Feige 


34 




JDbar 


0orb a 


V 


JDbar 


V 




-2450000 




(km s _1 ) 


-2450000 


(kms~ 




-127.199 


0.027 


-282 ±27 


-127.184 


43 ± 10 


44.040 


0.220 


-91 ± 11 


44.058 


40 ± 


4 


44.150 


0.401 


194 ±12 


44.161 


20 ± 


6 


46.131 


0.678 


74 ± 19 


46.156 


20 ± 


6 


46.146 


0.704 


43 ±10 








47.097 


0.275 


-19 ±19 


47.110 


13 ± 


4 


103.989 


0.362 


138 ± 7 


103.973 


7± 


4 


104.014 


0.404 


172 ± 7 









a Using the ephemeris of Lorimer et al. ( 1995 ) 
JD bar 2449220.947499(1), P orb = 0.604672713(5) d 



To 



but those taken with the 300 line mm -1 grating. For 
Feige 34, we found that the velocity derived from Ha 
was consistent with the results from the other spectra, 
while H/3, H7 and H<5 were off by 60 to lOOkms -1 . 
Therefore, we used the velocity from Ha alone for 
J1012±5307 as well. 

The velocities were corrected for possible small 
wavelength calibration errors using the O I A5577 
night sky emission line. The resulting velocities, cor- 
rected to the solar system barycenter, are listed in 
Table |^, and shown in Fig. ||. 

The radial velocities were fitted to a circular or- 
bit, with the period and time of ascending node fixed 
to the radio timing values of Lorimer et al. ( |1995| ). 
We find a radial-velocity amplitude ifwD = 265 ± 
9kms _1 and a systemic velocity 7 = —38 ± 6kms~ 1 
(Xicd = 0-97 for 6 degrees of freedom). If we ignore 
the uncertain velocity from the low-resolution spec- 
trum, K WB = 269±llkms^ 1 and 7 = -40±7kms" 1 
(X r 2 cd = 1-06 for 5 d.o.f.). 

The velocities determined for Feige 34 are, on av- 
erage, more or less consiste nt wit h the determinations 
by Greenstein & Trimble ( 1967 , 17 ± 13kms _1 ) and 
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Fig. 2. — Radial-velocity measurements for the com- 
panion of PSR J1012+5307. Overdrawn is the best- 
fit circular orbit, with the time of ascending node and 
period fixed to the values determined from the radio 
timing solution (Lorimer et al. 1995). The radial- 



velocity curve for the radio pulsar is shown as well 
(dotted curve). 
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Fig. 3.— The masses of PSR J1012+5307 and its 
companion. Shown are the constraints from the pul- 
sar mass function (thin, drawn curve for i = 90° ; dot- 
ted curves for inclinations as indicated), the ratio of 
the radial- velocity amplitudes (thick, slanted curve), 
and the white-dwarf surface gravity (thick, horizon- 
tal curve). For the latter two, the 95% uncertainty 
regions are indicated by thick, dashed curves. 



Thejll, MacDonald, & Saffer ( |1995| ; 3±4kms" 1 ). The 
individual velocities deviate significantly from the av- 
erage, however, indicating that some systematic ef- 
fects may be present. To estimate their effect on 
the orbital solution, we corrected the velocities for 
J1012+5307 with those for Feige 34 (corrected for 
an assumed velocity for Feige 34 of 5kms _1 ), and 
refitted the orbit. We find K W d = 290 ± Okms" 1 
and 7 = —62 ± 9kms _1 , with Xrcd = 0-73 (ignoring 
the low-resolution point: 297 ± 9, —66 ± 9 and 0.61, 
respectively). Here, the uncertainty in 7 includes 
a 5kms _1 uncertainty in the velocity of Feige 34. 
The Xrcd i s lower than for the uncorrected veloci- 
ties partly because the errors on the velocities now 
include the uncertainty in the velocities derived for 
Feige 34. Even if we neglect the latter, however, \rcd 
is still slightly smaller. At present, it does not seem 
possible to decide which method provides the correct 
answer, and therefore we will, conservatively, adopt 
if wd = 280 ± 15 km s" 1 and 7 = -50 ± 15 km s" 1 . 

From the radial-velocity amplitudes of the white 
dwarf and the pulsar (20.9774(2) kms -1 ; Lorimer et 



al. 1995 ), it follows that the mass ratio Mpsr/Mwd = 
13.3 ± 0.7. In Fig. ||, we show this constraint, as 
well as the constraints derived from the pulsar mass 
function (5.78361(10) 10~ 4 M Q ) and from the white 
dwarf mass. From this figure, it follows that with 
95% confidence the mass of the neutron star is within 
1.5 < M PSR /M Q < 3.2. 

5. Discussion and conclusions 

At present, the neutron-star mass determination 
is not accurate enough to constrain the equation of 
state, or to test the prediction that a large amount of 
mass has been accreted. To improve the accuracy, the 
uncertainties in both the mass ratio and the white- 
dwarf mass need to be reduced. The former will be 
relatively easy (and may well be interesting on its 
own; see Fig. ^), and we are planning more precise 
radial-velocity measurements for this purpose. 

It will be less straightforward to improve the white- 
dwarf mass determination. Although better models 
are becoming available (Hansen & Phinney, in prepa- 
ration), there are potential problems related to the 
uncertainty in the thickness of the hydrogen layer and 
the unknown abundance of helium in the atmosphere. 
Helium can be dredged up by convection if the hydro- 
gen layer is thin enough, and, if present, would mimic 
an increased surface gravity (Bergeron, Wesemael, & 
Fontaine |1991| ). Indeed, Reid ( |1996| ) finds that his 
mass determinations based on the observations of the 
Einstein redshift are systematically lower than the 
spectroscopic estimates for temperatures < 12000 K, 
which could be attributed to the presence of helium. 
IR observations could potentiall y reso lve this issue 
(Bergeron, Saumon, & Wesemael 1995 ). 

If an accurate distance could be determined (via 
timing parallax or otherwise), it could be combined 
with the flux and temperature of the white dwarf to 
yield a precise radius. If different from the radius 
predicted from our log g (0.028 ±0.002 Rq), one could 
assume either that the problem lies with the mass- 
radius relation (and not with helium pollution), and 
infer a mass from the radius and logg, or that there 
is helium pollution, but that the mass-radius relation 
is fine, and use that to derive a mass from the radius. 

From Fig. ^, one sees that the inclination i J> 60° 
for MpsR <^ 2 M Q . This raises the possibility of 
a model-independent, accurate determination of the 
white-dwarf mass and the inclination via Shapiro de- 
lay. With our precise mass ratio, this would allow one 
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to determine the neutron-star mass accurately. 

Finally, there is a small but measurable amount of 
light that the white dwarf will radiate due to repro- 
cessing of the pulsar wind at its atmosphere. This 
may provide yet another constraint on the geometry 
of the system. Fortunately, PSR J1012+507 is bright 
in the radio and the companion relatively bright in the 
optical, so that all these observations are possible. 

We thank Brad Hansen, Yanqin Wu, Fernando 
Camilo and Andrew Lyne for useful discussions. M.H.v.K. 
acknowledges a NASA Hubble Fellowship and S.R.K. 
grants from NSF, NASA and the Packard Foundation. 
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